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We have investigated the properties of a model of ID anyons interacting through a (5-function 
repulsive potential. The structure of the quasi-periodic boundary conditions for the anyonic field 
operators and the many-anyon wavefunctions is clarified. The spectrum of the low-lying excitations 
including the particle-hole excitations is calculated for periodic and twisted boundary conditions. 
Using the ideas of the conformal field theory we obtain the large- distance asymptotics of the density 
and field correlation function at the critical temperature T = and at small finite temperatures. 
Our expression for the field correlation function extends the results in the literature obtained for 
harmonic quantum anyonic fluids. 



I. INTRODUCTION 



For hard-core particles moving in two spatial dimensions, one can unambiguously define the notion of braiding of 
the particle trajectories by introducing the winding number n that gives the number of times the trajectory of one 
particle encircles another particle. This fact makes it possible to consider "anyonic" particles with fractional exchange 
statistics [l|, 0|, for which the wavefunction acquires the non-trivial phase factor e^'^'^'*, where k is the "statistical 
parameter", whenever n changes by ±1. This situation can be contrasted with the case of three spatial dimensions 
where one can define only permutations (no braiding) of point-like particles leading to only integer statistics, i.e. 
K = 0, 1 for bosons and fermions, respectively. In physical terms, the anyons in two dimensions can be viewed as the 
charge-flux composites for which the statistical phase arises as the result of the Aharonov-Bohm interaction between 
the charge of one particle and the flux of the other Q . Experimentally, anyons can be realized as quasiparticles of the 
two-dimensional (2D) electron liquids in the Fractional Quantum Hall Effect (FQHE) [?|. Individual quasiparticles are 
localized and controlled by quantum antidots in the FQHE regime [5,], and the transport properties of multi-antidot 
systems should provide direct manifestations of their fractional exchange statistics @ . Dynamics of individual FQHE 
quasiparticles attracted considerable attention (see, e.g., 7, 8]) as a possible basis for realization of the topological 
quantum computation [9]. 

Both conceptually and in practice (e.g., in FQHE systems), the 2D anyons can be confined to move in one dimension. 
There are, however, the aspects of fractional statistics in one dimension that make its introduction more complicated 
than in two dimensions. One is that for strictly ID particles, a trajectory of one particle can not wind around another, 
making the sign of the exchange phase g±*^«sgn(xi-a:j )/2 fj^g^^ wavefunction should acquire when the particle with 
coordinate Xi moves past the one with Xj, undetermined. The sign of this phase depends on whether Xi rotates 
clockwise or counter-clockwise around Xj in the underlying 2D geometry, which also explains why the signs of the 
phase change at Xi = Xj are opposite for the two particles in the pair: rotation of one sense for increasing coordinate 
Xi implies the opposite rotation for increa sing Xj. This fact hindered the early attempts at direct introduction of 
the ID anyons as charge-flux composites [Ifl, UM- If implies that any description of the ID anyons requires an 
additional convention on the choice of the sign of the statistical phase for each pair of particles. As discussed in more 
details below, this choice can be arbitrary and affects the appropriate boundary conditions of the quantum-mechanical 
wavefunctions of the system of anyons. 

Another complicating aspect of the fractional statistics in ID is the interplay between the two types of statistics, 
exchange statistics discussed in the preceding paragraph and the exclusion statistics defined through the volume of 
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the phase-space occupied by one particle [1^ • The exclusion statistics provides effective description of the dynamic 
interaction of particles, while the exchange statistics is associated with the "real" non-thermodynamic statistical 
effects that continue to exist in the limit of hard-core particles with infinite repulsion. The model of ID anyons with 
(5-function interaction considered in this work contains both types of effects, and the interplay between them can be 
seen in Eq. (j25[) for the renormalization of the dynamic particle-particle interaction by the exchange statistics. The 
renormalized interaction determines the thermodynamics of the model and can be expressed in terms of the exclusion 
statistics. However, in the hard-core limit c — ^ oo, effective interaction constant is essentially independent of the 
exchange statistics, and the main features of the thermodynamics of the model coincide with that of free fermions. 
There are still anyonic effects (e.g., the shift of the quasiparticle momenta by the parameter of the exchange statistics 
k) in this limit. 

The purpose of this work is to provide a systematic description of the ground state, low-lying excitations, and 
the asymptotics of the correlation function of the gas of ID anyons with the (5-function repulsion. In the form used 
below, the model was introduced by Kundu jlS] , who also provided the Bethe Ansatz solution. It was further analyzed 
recently by Batchelor et al. ffl [isL [l6l | . The model is an anyonic extension of the Bose gas with (5-function interaction 
solved by Lieb and Liniger [17|. The anyon gas reduces to the Bose gas in the limit when the statistical parameter k 
vanishes. In general, the Bethe equations for anyons are equivalent to the Bethe equations for bosons with two effects 
of the statistics k: renormalized coupling constant, and a twist in the boundary conditions. This makes it possible 
for us to use some of the known results for the Bose gas in the discussion of anyons. (Detailed description of the 
Bose gas with (5-function interaction including the correlation functions can be found in [l8|). The main results of our 
work are the formulae ([77)) . ([5^ for the density correlation function and ([57|) . ([5^ for the field correlation function. 
Expressions for the field correlators extend the results of Calabrese and Mintchev 19j| obtained in the harmonic fluid 
approach by including the higher-order terms that correspond to particle-hole excitations. Also, the conformal field 
theory approach we use provides immediate generalization of the zero-temperature correlators to finite temperatures. 

There are other ID models of anyons in the literature. Liguori, Mintchev and Pilo |2C|] investigated the momentum 
distribution of a more general gas of free anyons and predicted anyon condensation in a certain range of the statistical 
parameter. Ilieva and Thirring \2l\ studied the Hilbert space structure of the anyonic field, and showed that for a 
fixed statistical parameter it can be represented as an orthogonal sum of sectors with different numbers of particles. 
The Hilbert space of our model has the same structure. 

The paper is organized as follows. Sections [III introduces the field theoretical model of the ID gas of anyons with 
periodic and twisted boundary conditions. Equivalent quantum-mechanical problem is formulated in Section IIIII 
In Section IIV[ we discuss the properties of the ground-state, and in Section |V] calculate the finite size corrections 
for the ground-state and properties of the low-lying excitations. In Section IVIl using the conformal field theory 
approach we find the large-distance asymptotics for the zero-temperature density and field correlation functions and 
correlators at small finite temperatures. Appendix [X] presents the discussion of the boundary conditions for many- 
anyonic wavefunctions used in Section IIIII Appendix |B] describes the calculation of the energy and momentum of 
particle-hole excitations for periodic and twisted boundary conditions. 

II. THE LIEB-LINIGER GAS OF ANYONS 

We consider a gas of anyons with (5-function interaction in one dimension characterized by the Hamiltonian 

H= dx {[d.,^\{x)][d^<fAix)]+c^\{x)<f\{x)<fAix)<fA{x)}, (1) 
Jq 

where c > is the coupling constant and L the length of the system. The anyonic fields obey the equal-time 
commutation relations 

*A(a;i)*Ii(a;2) = e-'^'''^^'-^-^^\{x2)^A{xi) + 6{x, - x^) , (2) 

^\{x^)^\{x-,) = ^)*t^(2;2)vl/t^(xi) , (3) 

^ a{x^)^ a{x-,) = e*^''^(^^-^^)vl,^(a;2)*A(xi) , (4) 



where 



1 when x\> xi, 

e{x\ — x-i) = — 1 when x\ < xi, (5) 
when x\ ~ X2 ■ 
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In the original work introducing this model, the anyonic fields were realized in terms of the bosonic fields 
where 

p{x) = ^\{x)^A{x)^'fUx)'fB{x). (7) 

Due to the fact that at coinciding points e(0) = 0, the commutation relations (|2|3|4p are indeed bosonic. An alternative 
realization in terms of the fermionic fields was proposed in [2^ . However, in this case, the interaction term in the 
Hamiltonian ((T|) vanishes, since ^'^(a;) = [^'^(x)]^ = in coinciding points (see also the discussion in [l^). One 
implication of this difference is that in comparison to the bosonic representation ([6|) , similar fermionic representation 
with appropriate modification of the statistical parameter, effectively makes it possible to describe only the infinite 
repulsion limit c — > oo. 

Characteristics of the anyonic gas ([1]) depend on the boundary conditions imposed on the system &t x — Q ~ L. 
In this work, we use two different quasiperiodic boundary conditions which impose periodicity either directly on the 
anyonic or on the bosonic fields. Equations ([6]) imply that the periodic boundary condition for anyons correspond to 
twisted boundary conditions for bosons and viceversa. In terms of the anyonic fields, the boundary condition we use 
are, 

periodic BC: *^(0) = *^(L) , (8) 

and 

twisted BC: ^^^(O) = *t^(L)e"""(^"i' , (9) 

where N is the number of particle in the system. One can see directly from Eq. ^ that the external phase shift 
'kk{N — 1) introduced into the conditions (HI), ensures the periodicity of the bosonic fields. As will be shown in more 
details below, this means that this phase removes the anyonic shift of the quasiparticlc momenta. Below, we use the 
common notation for the two types of boundary conditions: 

■^\[0)^m\{L)e''^f^<^^^\ ^ = 0,1. (10) 

An important difference of the anyons with fractional exchange statistics from the integer-statistics particles is that the 
boundary conditions (jlOp for the fields do not translate directly into the same boundary conditions for the quantum- 
mechanical wavefunctions of the iV-anyon system p], which have more complicated structure (|23p derived in the 
Appendix VK\ 

The corresponding equation of motion —idt^A{x, t) — a{x, t)] for the boundary conditions (fTU]) is the nonlinear 
Schrodinger equation 

idt^Aix, t) = d^^A{x, t) + 2c^\{x, t)^\{x, t) . (11) 
The number of particle operator Q and the momentum operator P are defined as 

L 

dx -^^ix)-^ a{x) , (12) 



P 



I 

2 Jo 



L 



dx [^^Jx)d,^Aix) - [a,*^(x)]*^(x)) . (13) 

Both of them are hermitian operators which commute with the Hamiltonian 

[H,P] = [H,Q]=0. (14) 

If we define the Fock vacuum as 

^a{x)\0) , x£[0,L], (15) 
the A''-particle eigenstate of the Hamiltonian (and also of P and Q) can be then written as 

= / d^x e-'^XNixi, • • • , XN)'ifiixi) ■ ■ ■ 'i'i{xN)\0) , (16) 
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where the many-body wavefunction obeys 

XAr(a;i, • ■ • ,Xi,Xi+i, ■ ■ ■ ,xn) = e~''"^'''^'''~'''+^'> xn{xi, ■ ■ ■ ,Xi+i,Xi, ■ ■ ■ ,xn) ■ (17) 

This can be seen directly by using the exchange relation of the field operators '^\{xi)'^\{xi^i) — 
g'''^t^<^i^i-^i+i)'^'<^(^Xi+i)'^\^(xi) and interchanging the name of the integration variables Xi, Xi+i. Iterating the ex- 
changes several times we obtain 

III. THE EQUIVALENT QUANTUM MECHANICAL PROBLEM 

In [ll,[ii, it was shown that the eigenvalue problem (for periodic boundary conditions) 

H\i;)N = EN\^P)N, P\^P)n^PnH)n, (19) 
can be reduced to the quantum-mechanical problem 

H\ij)N ^ EnWn, V\^P)n^Pn\^}n, (20) 

where 

^ / fl2 \ 

^^^ = E +2c E ^(^.-^fc)' (21) 

j=l \ J / i<j<k<N 

(22) 

These considerations also hold for twisted and all cyclic boundary conditions for field operators. The boundary 
conditions for the quantum-mechanical wavefunctions of N anyons are (see 6] and Appendix [A|) 

Xjv(0, X2, • • ■ , xn) = e^-'^^^^-i) xn{L,X2,-- - , ^n) , 

Xjv(xi,0,--- ,a;jv) = e-'^^^e'^^'^^^-'^Nixi, L, ■ ■ ■ , xn) , 

■■ (23) 
Xn{xuX2,--- ,0) = e-''^'''''>-^e'-^-^^-'^XN{xi,X2--- ,L), 

where, as defined above, /? = 0, 1 for periodic and twisted boundary conditions pH)) . 

Using the Coordinate Bethe Ansatz [l^, 0, [3| we can obtain the eigenfunctions of the Hamiltonian pT|) as 

^ . E,<.^(^. J2i-ir^e^^"-^^-'-^ Y[[Xr,-Xv,-^c'eix,-x,)], (24) 



^^^!a->J(A,-Afc)2+c'2]- 
where [V] is the signature of the permutation and 

c' = 7^-^ (25) 

cos(7rK/2) ^ ' 

is the coupling constant renormalized by the exchange statistics. The eigenvalues of the Hamiltonian and momentum 
operators are = X^jLi Aj and = X^jLi Aji respectively. For the boundary conditions (j23p we obtain the Bethe 
equations 

N 



- n (28) 

k=i,k^j 



The Bethe equations (p6)) are similar to those obtained by Lieb and Liniger for the Bose gas with repulsive (5-function 
interaction. In our case, however, the effective coupling constant c' ([25|) can take negative values. While it can be 
shown (see, e.g., |18| ) that the Bethe roots are real for c' > , the roots can become complex for c' < 0, and one 
gets bound states [2^. In this work, we will consider only the case c' > 0. 
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IV. PROPERTIES OF THE GROUND STATE 



Bethe equations (|26)) can also be written as 

N 



AjL + ^6l(Aj - Afc) ==27mj+7rK(l-/3)(A^-l), j = l,---,N, (27) 



fc=i 

where 



and rij are integers when N is odd and half-integers when N is even. 

A. Twisted Boundary Conditions 

In this case {f3 = 1), the Bethe equations are similar to those for the Bose gas with periodic boundary conditions 
[13, [3 with c' as a coupling constant. The ground state is characterized by the set of integers (half-integers) 
rij — j — {N + l)/2, so the Bethe equations take the form 

XfL + f2e{\f -\^) = 2n f J , j^l,---,N. (29) 
fc=i ^ ^ 

From now on the superscript B will mean that the variables and physical quantities are the same as the ones for the 
Bose gas with periodic boundary conditions and coupling constant c'. In the thermodynamic limit N, L oo, D = 
N/ L =const, the Bethe roots become dense and fill the symmetric interval [— g, q\ . The density of roots in this interval 
obeys the Lieb-Liniger integral equation 

P(A) - ^ /" K{\ fi)p{f,) dfi=^, (30) 

where K{X,^) — 9'{X — = 2c'/(c'^ -I- (A — /i)^). The Fermi momentum q can be obtained from the Lieb-Liniger 
integral equation and the particle density is 

D = Y= j pW'^^- (31) 
Finally, the energy and the momentum of the ground state are 

r X^p{X) dX, Pq ^0. (32) 

J-q 

B. Periodic Boundary Conditions 

This is the case treated in [13, [H, [3 ■ The Bethe equations (|27)) in this case (/3 — 0) are similar to those for the 
Bose gas with twisted boundary conditions: 

N 

XjL + Y^e{Xj - Xk) = 2Trnj +ttk{N -I) , j = l,---,iV. (33) 

k=l 

Introducing the notation {[...]} such that 

{[a;]} =7, ifa; = 27rx integer + 27r7 , 7 £ [0, 1) , (34) 
we can describe the ground state by the following set of the Bethe equations: 

N 

A,L + ^0(A, -Afe)-2^(j-^^J +27rJ, j^l,---,N, (35) 
fc=i ^ ^ 
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where 5 = {['kk,{N — 1)]}. Comparison of Eqs. ([35]) and (|29p shows that we have the following connection between 
the Bethe roots for periodic and twisted boundary conditions: 

Xj =\f + 2t:5/L. (36) 

This relation is exact and holds also for the excited states if the (half)integers in the Bethe equations are the same. 
In the periodic case, the ground state is shifted by 2-k5/L, so that the Bethe roots are now distributed in the interval 
[— <7 + 2715 / L, q + 27r(5/L], and momentum of the ground state Pq in general does not vanish: 

N N 

Pa=^\,=Y^{\f + 2ti6/L)^2tiD5. (37) 

i=l 4=1 

The ground-state energy is: 

E„.t^^.tii^ff.'-^,^2^).ES.^, (38, 

i=l i=l ^ ^ 

where we have used that the total momentum in the case of twisted boundary conditions is zero and Eq in the 
thermodynamic limit is given by Eq. ([32|l . 



V. FINITE SIZE CORRECTIONS 



In this section, we are going to calculate the finite size corrections for the energy of the ground state and charac- 
teristics of the low-lying excitations. Based on the results of this section, we will be able to find the large-distance 
asymptotics of the correlations functions using conformal field theory. A chemical potential h is added to the Hamil- 
tonian ([1]) throughout this section, so that the total Hamiltonian is 

H„ = / dx{[5,vI;t^(x)][a,*^(a;)]+c*^(x)M'^(x)*^(a;)vI'^(x)-ft^*^(x)vI/^(x)}. (39) 
Jo 



A. Finite Size Corrections for The Ground State Energy 



As we have seen in the previous section, the ground state of the gas of anyons with twisted boundary conditions 
(/3 = 1) is characterized by the same set of Bethe equations as the Bose gas with coup ling constant c' and periodic 
boundary conditions. So in this case we can use the results for the Bose gas [Tsl. [23. [25l. l26l [27l [28j : 



= Lf £o(A)p(A) dA - ^ + O (^-1) , 



(40) 



-a 

where £o(A) = — h and vp is the Fermi velocity for the Bose gas with coupling constant c'. In the case of periodic 
boundary conditions (/3 = 0), Eq. then gives: 

Eo - L £ e„(A)p(A) dA - ^ + « + (i,) . (41) 



B. Low-Lying Excitations 



In our discussion of the low-lying excitations, we consider several different types of excitation processes: 

• Addition of a finite number AiV of particles into the ground state of the system. 

• Backscattering: all integers nj in the set {rij} characterizing the ground-state distribution are shifted by an 
integer d. 

• Particle-hole excitations: the integer rij that characterizes the particle at the Fermi surface is modified from its 
value in the ground state distribution by TV"*" for the particle with momentum q, (or q-\-2T:S/ L, depending on the 
boundary conditions) or by N~ at the opposite point of the Fermi surface with momentum —q, {—q + 2715 / L). 
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The central feature of the gas of anyons is that the boundary conditions for the field operators and the wavefunctions 
depend on the number of particles in the system. This means that any modification of the number of particles in the 
system changes the Bethe equations and, as a result, the quasiparticle momenta given by the Bethe roots. If we add 
one particle to the system of N particles, the boundary conditions are: 

XN+i{xi,0, - ■ ■ ,XN,XN+i) = e"'^'"'e*'''^'''-^"^^XAr+i(a;i, L, • • • ,xn,xn+i), 



(42) 



XN+iixi,X2,--- ,XN,0) = e-'^^^''e'^^-^^-'hN+i{xi,X2--- ,XN,L), 



and the Bethe equations become 

N+l , . . . , s 

giA.L ^ g*77«7Vg-i7r/3K(JV-l) TT / Aj - Afc + \ 

k=l,k^j V J / 

The ground states for N and + 1 particles are characterized by the Bethe roots satisfying different equations: 





N 

f^0(A, 






k=l 






N+l 




A,L4 


- E 

fc=i 


-Afc) = 27r(^j 



7V + 1 
2 

iV + 2 



27rw, j = l,---,7V, 

I-kJ, j = 1, • • • , at + 1 , (44) 



where 

w = 0, cj' = k/2, and cj = {[7rK(iV - 1)]}, J = {\t:kNW, 

for the twisted (/3 = 1) and periodic (/3 — 0) boundary conditions, respectively, and {[...]} is defined by Eq. (|34p . 
Comparing Eq. (gH) with Eq. we see that 

Aj = Af^r + , A, = Af + 2t:J IL , (45) 

where A^ are the Bethe roots characterizing the ground state of a gas of N bosons with periodic boundary conditions 
and coupling constant c'. 

1. Addition of One Particle to the System 

For excitations of this type we assume that both before and after the addition of a particle, the system is in the 
ground state. In order to calculate the energy and momentum of this excitation, we use Eq. (|45p which enables one 
to express energy and momentum through corrections to the same characteristics of excitations of the Bose gas. 

For the energy we get from Eq. P5|) : 



N+l N 

AEiAN = 1) = E ^o(^J-) - E 

^ AE^ {AN ^ 1) + {N + 1) i^^^ j , (46) 

where AE^JAN = 1) is the energy of the corresponding bosonic excitation. As known in the literature (see, e.g., 
[3 113, HE 113, nil) it is convenient to express this energy in terms of the "dressed charge" Z{X): 



(47) 
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where Z = Z{fi) = Z{—q), and Z{X) is defined as solution of the equation 



27r 

From (gS]) and (gT]) we obtain 



^W-i: I K{X,fi)Z{fi) dfi=l. (48) 

-9 



The momentum of the excitation is: 

N+l N 

^ \ . - (AT ^ ^\^ , 

L L ' 



AP(A7V = 1) = E - E = + - ' (50) 



where we again used the fact that for the ground state of bosons with periodic boundary conditions and any number 
of particles the total momentum is vanishing. 



2. Backscattering 

The uniform shift of the ground-state distribution in a backscattering process can be understood as a jump of some 
number d of particles between the opposite boundaries of the Fermi surface. The Bethe equations relevant for this 
process (in the case of N and N + l particles in the ground state) take the form: 

N / s 

k=l ^ ^ 

^'^+E^(^'-^fe) = 2^ J'-^— )+27rd + 2W, ,7V + 1. (51) 

fe=i ^ ^ 

Again, comparison with Eq. (p9|) shows that 

= + + d)/L , A^'' - Af^+i + 2niuj' + d)/L , (52) 
and the ground states are characterized by Eq. (|45|) . Using Eqs. (|45p and ((52| we get the excitation energy: 

N particles: I^E{d) = j^ie^iXj) - So(A,)) = iV - , (53) 



iV + 1 particles: ^E{d) = ^ (^^(A^^) - ^^(A,)) = (A^ + l)i^^^^^^^±^ - (A^ + 1)^^^ 



L2 ^ ' ^ L2 



(54) 



This result can be rewritten using the relation = 27^0 /vp (see [T^, Chap. 1.9) obtaining 

iV particles: /^E{d) = ^^^^(d + cj)^ - ^^Z^w^ 

iV + 1 particles: A£:(d) = i^Z^ d + Jf -Z^oj'^ + ^ ^ - ^ . 55) 

L L L'^ 

The momentum of the backscattering excitation is simply 

AP(d) = N{2nd/L) , (56) 
the expression that is valid for any number of particles N. 
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3. Particle-Hole Excitations at the Fermi Surface 

In this case, the excitations we consider consist in changing the maximal (minimal) nj in the ground state by 
For N particles and "excitation magnitude" the Bethe equations are 



XjL + J20i\j-\k) = 27t(j-^^) +2TTU;, j = 1, ■ • • , iV - 1 , 

k=l ^ ^ 

Aa,L + ^6I(Aa, - Afc) = 2tt (n - ^^^^j + 2nLj + 2nN+ . (57) 

k=l ^ ^ 

Prom ([57]) we see that the momentum of the excitation N'^ is AP(iV+) — 2ttN~^ /L and, similarly, for the excitation 
N~ the momentum is AP{N~) = —2ttN~/L. These excitation can be considered as a special case of the general 
particle-hole excitations, and we can use the results of Appendix [Bl for them. Using (jB7|) we see that the excitation 
energy and momentum 

Ai?(iV±) = ^iV± + O (^-1) , AP(iV±) = ±^7V±, (58) 
coincide with those for the similar excitations of the Bose gas (see Appendix 1.4 of p^): 

A£;^(7V±) = ^7V±, AP^(7V±) = ±^7V±. (59) 
For + 1 particles, the energy and momentum of the excitations are given by the same expressions as in (j58p . 



VI. LARGE-DISTANCE ASYMPTOTICS OF CORRELATIONS FUNCTIONS 



In this section we calculate the asymptotics of the correlation functions. We will consider the case of twisted 
boundary conditions (/3 = 1), or the periodic boundary conditions (/3 = 0) when k is a integer multiple of 2/{N — 1), 
so that the shift in (|36p vanishes, ^ = 0, and the two boundary conditions are equivalent - see (fTO|) . The main 
feature of this case that is important for the direct applicability of the conformal field theory approach is that the 
momentum of the ground state p7p of the gas of anyons is zero for these boundary conditions. For general gapless 
1+1-dimensional systems, T = is a critical point making the correlation functions decay as a power of distance at 
T — but exponentially at T > 0. As we have seen in the previous section, the Lieb-Liniger anyonic gas is gapless 
and the excitation spectrum has a linear dispersion law in the vicinity of the Fermi level. These features support the 
expectation that the critical behavior of the anyon system is described by conformal field theory (CFT). 

CFT is a vast subject and we refer the reader to [1^ H^, [3l|, and Chap. XVIII of for more information. 
A conformal theory is characterized by the central charge c (not to be confused with the coupling constant in ^) 
of the underlying Virasoro algebra, and conformal invariance constrains the critical behavior of the systems under 
consideration. The critical exponents (the powers that characterizes the algebraic decay at T = 0) are related to 
the conformal dimensions of the operators within the CFT, so to obtain the complete information about the critical 
behavior of the system we need to calculate the central charge and the conformal dimensions of the primary fields. 



A. Central Charge 

In order to find the central charge we use the fact that for unitary conformal theories it can be found from the 
finite-size corrections, specifically the coefficient of the 1/L term in the expansion of the ground state energy for 
L-^oo [Hill: 

i.^X..-^c + o(i). (60) 

Comparing this relation to Eq. (I40p valid for the boundary conditions we are assuming in this Section, we see that the 
central charge c = 1. The fact that the central charge c = 1 means that the critical exponents can depend continuously 
on the parameters of the model [2^, [s^ [s^ . 
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B. Conformal Dimensions from Finite Size Effects 

Following the original idea of Cardy |32] subsequently developed in [23, H^l , we obtain below the conformal 
dimensions of the conformal fields in the theory from the spectrum of the low-lying excitations described in the 
previous Section. The local fields of the model can be represented as a combination of conformal fields 

(61) 

Q 

where A(Q) are some coefficients and z = ix + vpT, with vp the Fermi velocity and t the Euclidean time. The 
conformal fields are related to excitations with quantum numbers Q = {AN,N^,d}, where AN represents the 
number of particles created by the field cj), and all the fields in the expansion ((6T|) should have the same AA^. The 
quantum number d gives the number of particles backscattered across the Fermi "sphere" , and characterizes the 
change of the maximal or minimal rij in the Bethe equations from its values in the ground state. While A7V has to 
be the same for all the terms in the expansion, d and can be different. 

For two conformal fields, (j)Q and (j)Q' , with the same conformal dimensions denoted A^, their correlation function 
is given by 

{(j)Q{zi,Zi)(j)Q,{z2,Z2)) ^ ^^A^ N2A- • (62) 

(Zl - Z2y^ (Zi - Z2y^ 
Under a conformal transformation z = z(w), z — z{w), it transforms like 

^ f dz f dz 

{(I)q{wi,Wi)(I)q,{w2,W2)) (^a^j ('?!'q(^:i(i«i),^i(wi))</'Q'(2:2(w2),^2(w2))) ■ (63) 

Using the expansion (|6ip . the fact that the two conformal fields with different conformal dimensions are orthogonal 
(their correlation function is zero), and (|62p we then have: 

(^(zi,z-i)</>(z2,.-2)) = (64) 

Q (21-Z2) «(2l-^2) « 

which is valid in the whole complex plane without the origin (zi 7^ 22)- Conformal mapping of this plane to a cylinder 
(periodic strip) with the help of transformation 

^^g2Wi^ w = ix + VFT with 0<2;<L, (65) 

applied to ([63]) gives 

(0(zz;i,«;i)0(^2,«)2)) -E^(Q) f • M wn V'''' f-TTT^^^^^^TTTiV^'' ' (66) 

with the asymptotics 

((/.(w;i,M)i)0(«;2,«;2)) -^e-^(^«+^Q)("i-"^)-'T(Aj-AQ)(.i-x2)^ (67) 
Q 

Comparison with the spectral decomposition of the correlation function in the periodic strip (ti > T2) 

(0(w;i,u;i)</.(u;2,t<)2))L = ^|(O|0(O,O)|g)pe-(^«-^°)(^i-^^)-'(^«-^«)("i-"^), (68) 

Q 

where |0) is the ground state and Eq , Pq are the energy and momentum of the ground state, leads to 

Eq^E, = ^(A+ + Aq) , - P„ = ^(A+ - Aq) , (69) 

assuming that both the energy and momentum gaps are of order 0{\/L). However, as we have seen in Sect. |Vl 
for some of the excitations considered (addition of a particle in the system, AiV = 1, backscattering processes 
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characterized by d, and particle- hole excitations at the Fermi surface characterized by N^), the momentum gap is 
macroscopic. For example, ii Q = {AN — 0,d ^ 0, = 0}, the momentum gap is 2kpd, kp = t tD, and for 
Q = {AN = l,d ^ 0, N^ = 0} the momentum gap is nkpK + ttk/L. For these excitations, following [23, 0, H^j the 
coefficients A{Q) will depend on x as 

A{Q) = A(g)e'P«^ , (70) 

where pg is the macroscopic part of the momentum gap Pq — Pq. From (j64p and (j70p we obtain the generic formula 
for the asymptotics of correlations functions at T = 

t)^(0, 0)) = ^ ' (71) 

where Ag can be found from ([M)) and the leading term corresponds to the smallest Ag. 

We also can find the low-temperature asymptotics of the correlation functions if we use instead of the conformal 
mapping (|65p. the mapping 

z = £2^^'"/^^ , z^x- ivpT , (72) 

which differ from (|65p by interchanging the space and time variables. The computations are similar those described 
above for the correlation functions in a finite box, and the final result is 

\smh[TTi [x — tvpT)/vp\ J \sinh[ni [x + iVFT)/vp\ J 
This result is valid only at temperatures close to zero. 

C. Density Correlation Function 

In the case of the density correlation function, (j(x, t)j(0, 0)), where j{x) = '^\{x)'^ a{x) , we have AN = so the 
most general excitation is constructed by backscattering d particles and creating a particle-hole pair at the Fermi 
surface characterized by N^ . Making use of (|55l56l58p , we obtain for the energy and momentum gap of the excitation 
characterized by Q {AiV ^Q,d,N^}: 

PN±,d-Po = 2kpd+'^{N+ -N-), (74) 

EN±,d -Eo = ^ [{Zdf + N+ + N-]. (75) 

Here we have taken into account only the terms of order 1 and 0{1/L). Equation gives the conformal dimensions 

2A^ = 2N^ + {Zdf , (76) 

and from the general formula ((7T|) 

J2ixkFd 

0-(x,t)j(o,o)}-(j(o,o)}^= ^(Q) iX-- (77) 

Q={N±,d} (tx + vpr) Q{-ix + vpT) Q 

Defining 6 = 2Z^ = AttD/vp, where Z = Z{—q) = Z{q), and Z{X) given by the integral equation (|48l) . the leading 
terms are 



{j{x, t)m 0)) - m or - / + — ^ + . (78) 

(iX + VpTy (—IX + VpT)^ \lX + VpT\" 

For equal times, Eq. ([77)) takes the form 



^2ixkFd 



0-(x,0)j(0,0))-(j(0,0))2= Y ^(0)7-^i^T51vTT51v-- (79) 
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The presence of the oscillatory terms in this expression can be explained by the following simple computation [27[: 

0-(x,o)j(o,o)) = ^(o|j(:^:,o)|g}(g|j(o,o)|o) = ^|(o|j(o,o)|g)|V(^«-^°)- = 

Q Q 

OO 

= e'2''^''"5]l(0|j(0,0)M,A^±)pe"^(^^-^"), (80) 

d=-oo Af± 

where in the second line, we broke the sum over Q into disjoint sums characterized by different macroscopic momenta. 
The second part of the sum gives the power-law decay for kp^ <^ x <^ L. The formulae (|77|) and ((79| are the same 
as in the case of a Bose gas with coupling constant c' = c/ cos{ttk/2) and periodic boundary conditions [13] - see 
Chap. XVII of fisl]. This situation is expected, since 

jix) = ¥^ix)-^A{x) - ^U^)^Bix) . (81) 
Finally, from (j73p . the finite temperature density correlation function is 

(,(.,.),(0,0)).= 5: BiQ)e-^^^^ ( T U -^^T^i'" V ^''^ 

J \smh\TT 1 [x — ivft)/vf\ J \smh\Trl [x + ivft)/vf\ J 

with Ag given by ([751) . 



D. Field-Field Correlator 



In contrast to the density correlators, for the field correlator (^yi(a;, i)5'^(0, 0)}, one has AA^ = 1. For the ground 
states with N and + 1 particles and the boundary conditions considered in this Section the Bethe equations are: 



N 



k=l 

N+1 



A,L+;^^^(A, -Afe) = 2^(j-^^] j = l,-.. ,iV+l. (83) 

k=l ^ ^ 

The shift ttk in the second equation implies that the anyonic wavefunctions for N and N + 1 particles live in two 
orthogonal sectors of the Hilbert space. The addition of one particle produces in this case a macroscopic change in 
the momentum, TrkpK + ttk/L, which gives rise to oscillations even in the dominant term of the field correlator. 

The most general excitation is obtained by an addition of one particle to the system, followed by the backscattering 
of d particles and creation of a particle-hole pair at the Fermi surface. Using the results (|49I50I55I56I58P with 
Lu = 0,Lu' — k/2, we obtain the following expressions for the energy and momentum gaps of an excitation with 
Q — {AN — l,d,N^} (retaining, as before, the terms of order 1 and 0{1/L)): 

Pn±,7^ 'Po^ 2fcF(d + k/2) +'^[id + k/2) + N+ ~ N'] , (84) 



2'kvf 
L 



^] +Z^{d + n/2)^ + N+ + N- 



(85) 



so the conformal dimensions are 



2A^ = 2N^ + [ ^ ±Z{d+ k/2] 



1 

2Z 



From Eq. (|7ip . the field correlator is 



={N±,d} (ix + vft) '^'^Q{-ix + VFT) 



(86) 



(87) 
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or in the equal-time case 

Q={N±.d} l-^l 

where 6 = 2Z^. A gain, we can heuristicaUy justify the presence of the osciUatory terms in the correlation function 
in the same way as for the density correlator, but for the field correlator, the complete set of states that is inserted 
between ^Pyi and ^'^ is from the sector with + 1 particles 



{^Aix, 0)*^(0, 0)) = ^(0|vI.^(.T, 0)\Q){Q\^\{0, 0)|0) = ^1(01*^(0, 0)|Q)| 



= J2 e'2fe^('^+t)-^|(0|*A(0,0)|d,iV±)|2e'^(^^"^"). (89) 

d=-oo N± 

In this case, the terms of the correlation function containing g«2/cF(d+K/2)2; fi^^Q^ ^^.q responsible for the oscillatory 
behavior at a; ^ L, exhibit dependence on the statistical parameter. 

Equation ([57)1 can be compared to the result of Calabrese and Mintchev 19], who calculated the field correlation 
function for anyonic gapless systems in the low- momentum regime using the harmonic fluid approach (38l.[39|. obtaining 

°° p-2i{d+^)kFX -2i{m+§)ire(x)/2 

K(x,0)vl.,(0,0)) , (90) 

where D is the density, bd unknown non- universal amplitudes, c{x) = L sin(7ra;/L), and X is a universal parameter that 
can be expressed in terms of the phenomenological velocity parameters vn,vj as K = \fvJJ%r^. For the Lieb-Liniger 
anyons, 

2-kD e , , 

91 

VF 2 

They have checked their results in the limit c — > oo, K = 1 against the exact results of Santachiara et al. [40*1 , who 
calculated the generalization of Lenard formula [4l| for anyonic statistics. We see that our conformal field theory 
approach agrees with the leading asymptotics produced by the harmonic liquid approximation but also gives the 
higher-order terms in the large-distance expansion. 

Using the conformal mapping (TfS]) that leads to general Eq. ((73|) , we find also the finite-temperature field correlator: 

J Xsmtiy-Kl {x — ivpT) I vp\ J ysmhy-Kl {x + ivft)Ivf\ J 

+ . . (^2) 

where Aq is given by 



VII. CONCLUSIONS 



The main result of our work is the calculation of the large-distance asymptotics of the correlation functions of the 
gas of ID anyons using the ideas of conformal field theory. This result requires conformal invariance close to the 
critical point T = 0, and the knowledge of the finite size corrections to the energy and momentum of the ground state 
of the gas due to low-lying excitations. In the analogous case of Bose gas with (5-function repulsive interaction, the 
conformal field theory predictions for the asymptotics of the correlators were checked against the exact results for these 
asymptotics obtained from the determinant representations and the differential equations for the correlation functions 
. It would be interesting to have similar exact results for the model studied in this paper which is a natural anyonic 
extension of the Bose gas. As a first step in this direction, Santachiara, Stauffcr and Cabra [13], already obtained for 
the one-particle reduced density matrix (field correlator) in the impenetrable limit a representation in terms of the 
determinant of a Toeplitz matrix of dimension (A^ — 1) x (A^ — 1) where A' is the number of particles. The exact results 
for the anyon correlation functions would also be needed to extend the correlators derived in this work for essentially 
one type of boundary conditions to more general quasiperiodic conditions. This problem seems particularly natural 
for anyons for which the effective boundary conditions for quasiparticle momenta change with the total number of 
particles in the system. 
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APPENDIX A: BOUNDARY CONDITIONS FOR THE MULTI-ANYON WAVEFUNCTIONS 

In this Appendix, we derive the exact form of the cycHc boundary conditions for the wavefunctions of the many- 
anyon system. Our treatment generahzes the approach of @ to the case of several penetrable particles. In physical 
terms, the situation we consider corresponds to anyons confined to move along a loop with, in general, an external 
phase shift (j) created, e.g., by a magnetic filed threading the loop. We start with the case of two particles and no 
external phase shift, (/> — 0. The Bethe-Anzatz wavefunction reduces in this case to the following form: In the 
region I (2:1 < X2) one has 

Xi{xux2) = -=__==={e'(-^^^+-^^^)(A2 - Ai - ic') + e*(-^^^+-^^^)(A2 - Ai + ic')} , (Al) 
V2[(A2 - Ai)-^ + c'^J 

and in the region II (xi > X2): 

Xii{xi,X2) = =={e-(-^^^+-^^^)(A2 - Al + zc') + e»(-iA.+..A,) (^^ _ ^ ^c')} . (A2) 

The general exchange symmetry of this wavefunction given by Eq. (|17p imply that for fractional k it can not satisfy the 
same boundary conditions in the two coordinates. As one can see by exchanging the coordinates, if the wavefunction 
is periodic in the first one, the boundary conditions in second one should have a twist, 

x{0,x)^x{L,x) ^ x{x,0)^x{x,L)e-^'^^, (A3) 



and viceversa. One consequence of this is that the exact form of the Bethe equations ([26|) depends on whether 
we impose periodic boundary conditions on one or the other coordinate. Indeed, if one requires periodicity in xi, 
x{0,X2) — x{L,X2), the Bethe equations are: 



whereas the periodicity in X2, x{xitO) — x{xi,L), results in the equations that differ by the sign of the statistics 
parameter k: 



n f^t^V («) 



Since the Bethe equations determine the spectrum of the quasiparticle momenta A^ through Eq. (|27p . the k shifts of 
different signs produce two physically different situations. 

The origin of this difference can be traced back to the fact that the fractional statistics requires braiding of particles, 
something that strictly speaking can not be done in one dimension. To define the braiding of ID particles one needs 
to first adopt a convention on how the particles pass each other at coinciding points, something that is done by 
choosing a specific sign of the exchange phase e~*'^'"^(^i~^2)/2 After that, one more choice that needs to be made is 
how the ID loop with anyons is imbedded into the underlying 2D anyonic system. In the case of two particles, this 
choice is reflected in the possibility of choosing different boundary conditions for two different anyonic coordinates 
and determines how the particle trajectories enclose each other as the particles move along the loop [6| . As reflected 
in Eq. (jA3p . periodicity in xi means that the trajectory of xi does not enclose the particle X2- This implies that xi 
is itself enclosed by the trajectory of X2, producing the twist in the boundary condition for X2 variable. The different 
choice of the boundary condition would mean that the ID loop in imbedded into the 2D system in such a way that the 
trajectory of xi encloses X2- This means that the wavefunction periodicity in both variables correspond to different 
but valid physical situations. 

The situation is somewhat more complicated for larger number of particles, as can be seen in the case of three 
particles. In the wavefunction (|24p . one needs to distinguish then six regions corresponding to the six permutation of 



the particles. The wavefunction ([24]) in these regions is: 
Region I (xi < X2 < X3) 



Region II {xi < x^ < X2) 



Region III {x^ < xi < X2) 

Xlllixi, X2,X3) 



Region IV {x^ < X2 < xi) 

XIV{XI,X2,X3) 



Region V {x2 < xi < Xz) 

Xv{xi,X2,x:i) 







-A2- 


- ic')(A3 - 


-Ai- 


- jc')(A2 - 


-Ai - 


- ic) 




_gi(£Ci Ai+a;2A3+a:3A2) _ 


A3 - 


ic')(A2 - 


Ai - 


ic')(A3 - 


Ai - 


ic') 




_|_gi(£ciA3+a;2Ai+i:3A2) _ 


Ai - 


*c')(A2 - 


A3- 


ic'){X, ~ 


A3- 


ic') 




_gj(£!;iA3+a:2A2+2;3Ai) _ 


A2- 


ic')iXi - 


A3- 


ic'){X2 - 


A3- 


ic') 




_|_gi(£CiA2+a:2A3+a;3Ai) j-^^ _ 


A3- 


ic'){Xi - 


A2- 


zc')(A3 - 


A2- 


tc') 




_gi(a;iA2+a;2Ai+a;3A3) j-^^ _ 


Ai- 


zc')(A3- 


A2- 


zc')(Ai - 


A2- 


J 




^ |g*(^l-''l+^2A2+2;3A3) j-^^ _ 


A2 -+ 


zc')(A3- 


Ai - 


*c')(A2 - 


Ai- 


ic') 




_ i(2;iAi+i:2A3+2:3A2)('\ _ 


A3 4 


ic')(A2 - 


Ai - 


ic')(A3 - 


Ai - 


ic') 




A3+2;2Ai+2:3A2) _ 


Ai -1- 


ic')(A2 - 


A3- 


^c')(Al - 


A3- 


ic') 




A3+2;2A2+2:3Ai) j-^^ _ 


A2 -1- 


ic')(Ai - 


A3- 


^c')(A2 - 


A3- 


ic') 




A2+a;2A3+2:3Ai) j-^^ _ 


A3 + 


ic')(Ai - 


A2- 


zc')(A3 - 


A2- 


ic') 




A2+2;2Ai+a;3A3) _ 


' Ai -f 


-»c')(A3- 


'A2- 


-ic')(Ai - 


-A2- 


- ic') 


Ae 


|gi(a;iAi+a;2A2+2:3>'3) 


-A2 


+ *c')(A3 


-Ai 


+ ic')(A2 


- Ai 


-ic 






- A3 ^ 


- ic')(A2 - 


- Ai 4 


-«c')(A3 - 


- Ai - 


- ic j 




A3+a;2Ai+2:3A2) 


-AH 


-*c')(A2- 


- A3 ^ 


-jc')(Ai- 


-A3- 


-zc') 




A3+a;2A2+i;3Ai) 


- A2 4 


- ic){Xi - 


- A3 4 


-*c')(A2- 


-A3- 


- ic ) 




A2+X2A3+2;3Ai) 


-A3^ 


-*c')(Ai - 


- A2 4 


-*c')(A3- 


-A2- 


-ic') 




_gi(xiA2+a;2Ai+a:3A3) 


-Ai- 


f«c')(A3- 


-A2- 


^^c')(Al - 


-A2- 


-ic') 


Ae 


=i|li |gj(a:iAi+i:2A2+2:3A3)^;^^ 


— A2 


+ *c')(A3 


-Ai 


+ IC')(X2 


- Ai 


4 ic 




_^i{xiXl+X2X3+X,,\2) ^^^^ _ 


- A3 -i 


- «c')(A2 - 


-AH 


-tc')(X3 - 


- Ai 4 


- ic) 




_|_gi(a;iA3+a;2Ai+a;3A2) _ 


-Ai4 


- «c')(A2 - 


-A3 4 


-*c')(Ai - 


- A3 4 


- ic) 




_gi(a;iA3+a;2A2+a:3Ai) _ 


- A2 4 


- «c')(Ai - 


- A3 -f 


-«c')(A2- 


- A3 4 


- ic) 




_|_gi(a;iA2+a;2A3+a;3Ai) _ 


-A3 4 


-»c')(Ai - 


- A2 4 


-ic')(A3- 


- A2 4 


-ic') 




_gi(a;iA2+2;2Ai+K3A3)^^^ _ 


-AH 


-ic')(A3- 


-A2H 


-«c')(Ai - 


- A2 H 


-ic') 


Ae- 


^ /g'(^l-^l+^2A2+2;3A3) _ 


A2- 


ic')(A3 - 


Ai- 


zc')(A2 - 


Ai + 


ic) 




_gi(a;iAi+K2A3+2;3A2) _ 


A3- 


*c')(A2 - 


Ai - 


ic')(A3 - 


Ai + 


tc') 




_|_gi(a;iA3+i:2Ai+2;3A2) _ 


Ai- 


*c')(A2 - 


A3- 


*c')(Ai - 


A3 4- 


tc') 




_gi(a;iA3+a:2A2+2;3Ai) j-^^ _ 


A2 - 


*c')(Ai - 


A3- 


*c')(A2 - 


A34 


tc') 




_^^i{xi\2+X2X3+X3\l) ^y^^ _ 


A3- 


*c')(Ai - 


A2- 


*c')(A3 - 


A2 4 


tc') 




_gi(a;iA2+a;2Ai+a;3A3) _ 


Ai- 


zc')(A3 - 


A2- 


*c')(Ai - 


A2 4 


tc') 



}■ 
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Region VI {x2 < X3 < xi) 

XVl{xi,X2,X3) 



I 


+X3A3)CX _ 


- A2 - 


-«c')(A3 - 


- Ai 


+ ic')(A2 - 


- Ai 


1 ■ 

+ ic 


i( X^ Ai A-X'> X'i A 


^^'\\2 — 


A3 - 


IC )(A2 - 


Ai 4 


ic )(A3 - 


Ai 4 


IC ) 


gi(a;iA3+a;2Ai-f 




Ai - 


ic')(A2 - 


A3 4 


*c')(Ai - 


A3 4 


ic') 


gi(a;iA3+a;2A2-f 




A2- 


ic')(Ai - 


A3-F 


*c')(A2- 


A3 4 


ic') 


gi(a;iA2+a;2A3-f 


^^^i)(Ai - 


A3~ 


ic')(Ai - 


A2 4 


zc')(A3- 


A2 4 


ic') 


_gi(2;iA2+2;2AH 


^^3A3)(^g - 


Ai - 


ic')(A3 - 


A2 4 


-*c')(Ai - 


A2 4 


- ic'\ 



(All) 



'}■ 



where 



A- 



6n,>J(A,-A.)2 4c'2] 



(A12) 



As discussed above for the two particles, the periodic boundary conditions can be imposed in principle on any of 
the wavefunction arguments. Requiring xi to be periodic, x(0, a;2,a;3) — x{L,X2,X3), gives 



Xl{0,X2,X3) 
Xll{0, X2,X3) 



XVliL, X2,X3) 
Xlv{L, X2,X3) 



for a;2 < a;3 , 
for X3 < a;2 . 



(A13) 
(AM) 



Except for the exchange-statistics phase factors, the wavefunctions in the six regions coincide with the wavefunctions 
of the Bose gas with the (5- function interaction of strength c' (|25p . Therefore, the Bethe equations we obtain are the 
same as in the bosonic case with the only difference coming from the statistical phase factors. Conditions (jA13p and 
(|A14[) produce six equations each, with only three of them being independent 



iL\2 



2i7TK 



2i7TK 



Ai 


-A2 


4 ic' 


Ai 


-A2 


— ic! 


A2 


-Ai 


4 id 


A2 


-Ai 


— ic! 


A3 


- Ai 


4 id 


A3 


-Ai 


— ic! 



Ai — A3 



Ai — A3 — ic 

A2 - A3 4 ic 



A2 - A3 - ic 
A3 - A2 4 ic 



A3 — A9 



(A15) 



These equations can be written in the compact form similar to Eq. (|26p : 

3 



iL At 



n 

k=l,k=^j 



A," — Afe 4 id 



A, 



Ai, - id 



If the periodic boundary conditions are imposed on the second variable, xi^ii Oi ^3) 

Xv{0,X2,X3) = xiiiL,X2,X3) , for xi < a;3 , 
Xvi{xi,0,X3) = xiii{xi, L,X3) , for xi > 0:3 , 

we obtain either from (|A17|) or (|A18p the following Bethe equations 

3 



xixi,L,X3), i.e.. 



n 



Aj - Afc + id 
Xj - Afc - ic' 



(A16) 



(A17) 
(A18) 



(A19) 



Finally, if we impose periodic boundary conditions on the third variable, xi^iT ^2,0) = xi^ij ^2, L), i.e., 

Xiiiixi,X2,0) = xi{xi,X2,L) , for xi < 2:2 , (A20) 



X/y(a;i,a;2,0) = xv{xi,X2,L) , for X2 < xi , 



the resulting Bethe equations are 



n 

k=l,k^j 



Xj - Afc 4 id 
Xj - Afc - id 



(A21) 



(A22) 
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The difference between the three forms of the Bethe equations (jA16p , (|A19|) , (|A22p means that the periodic boundary 
conditions imposed on one variable automatically require the twisted boundary conditions on the other variables if 
one wants to keep the same Bethe equations. Similarly to the case of two particles, this can also be seen directly from 
the anyonic exchange symmetry (|17p of the wavefunction. Suppose we set the periodic boundary conditions on the 
first variable: 

X{0,X2,X3) = x{L,X2,X3) . (A23) 

Exchanging then the first two variables on both sides of Eq. (|A23p with the help of Eq. (fT7|) . we get the twisted 
boundary conditions for the second variable: 

x(a:2,0,X3) = x(a:2,i,a:;3)e~2--. (A24) 

From (|A24p . using again (|17p we have 

Xix2,x3, 0) = x{x2,X3, L)e-4"'= , (A25) 

which are the twisted boundary conditions for the third variable which follow from the periodic conditions on the 
first. From any of the boundary conditions (jA23[l . (jA24|, (jA25P we obtain the Bethe equations (jA16| . 

Similarly, periodic boundary conditions on the second variable give the following boundary conditions for the 
three-anyon wavefunction: 

X{0,X2,X3) = x{L,X2,X3)e^"'" , 

x{xi,0,X3) = x{xi,L,X2), (A26) 
x(a;2,X3,0) = x{x2,X3,L)e~'^''"^ , 

and the Bethe equations (jA19p . The same can be done starting with periodicity in the third variable. As in the case 
of two particles, we sec that imposing periodic boundary conditions on the first and the last variables produces the 
Bethe equations, (|A16P and (jA22p . which differ only by the sign of the statistical parameter k. As discussed in detail 
for the two particles, this difference corresponds physically to different imbedding of the ID loop of anyons into the 
underlying 2D system. In the two situations, the number of particles enclosed by the trajectories of successive particles 
Xj, j = 1, 2, N, either increases from to — 1 or decreases from A^ — 1 to 0, as reflected in the corresponding 
boundary conditions of the multi-anyon wavefunction. In contrast to this, the requirement of periodicity of one of the 
"internal" variables (e.g., X2 in the case of three particles) produces the Bethe equations and boundary conditions, e.g. 
(|A19p and (|A26p . that do not have this interpretation. They describe the situations with appropriate non- vanishing 
external phase shift 7^ 0, which twists uniformly the boundary conditions of all the variables. In the main text of 
our paper, we use the periodic boundary conditions with respect to the first variable of the anyonic wavefunction or 
introduce the external twist (j) = tik{N — 1) which removes the anyonic shift of the quasiparticle momenta. As follows 
from the discussion in this Appendix, the boundary conditions for the wavefunction of N anyons are given in these 
two situations by Eqs. ((23|l . 

APPENDIX B: PARTICLE-HOLE EXCITATION 

In this appendix we find the energy and momentum of particle-hole excitations of the gas of anyons. As discussed 
in the main text, for twisted boundary conditions (/? = 1), the ground state of anyons is equivalent to that of the Bose 
gas with periodic boundary conditions and coupling constant c' , so the excitation energy and momentum coincide in 
this case with those known for the Bose gas (see Chap. 1.4 of fl8|). For periodic boundary conditions (/3 = 0), the 
Bethe equations are the same as for the Bose gas with the boundary conditions twisted by the phase shift 27r5, where 
5 = {[^^{N — 1)]}. In the case of one hole with momentum A/i and one particle with momentum Ap the equations for 
the ground state and the excited state are: 



JV 



Ground State, PBC : XjL + ^e{Xj - \k) = 2ti [ j —]+2ti5, j^l,---,N, (Bl) 



1 



N 



Excited State, PBC : XjL+^e{Xj-Xk)+e{Xj-Xp)-9CXj-Xh) ^ 2n (j - ^^^^]+2TrS , j = l,---,N. (B2) 



fc=i 
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Comparing the equations for a particle-hole excitation in the case of twisted boundary conditions 

Ground State, TBC : Af L + ^ 6l(Af - Af ) 27r U - V j = l,---,iV, (B3) 

k=l ^ ^ 



Excited State, TBC : Af L + ^ 0(Af - Af ) + 0(Af - A^) -0(Af - Af ) = 2^ (^j - , j = l,-..,iV, (B4) 

with (jBl[) and (jB2p . we find the following relations 

Aj = Af + 27r(5/L, A^ = Af + 27rVi, (j = l,---,iV) (B5) 



Ap = A^ + 2ttS/L, Xh =Xh + 2ttS/L . (B6) 
The energy and momentum of this excited state with respect to the ground state is (eq (A) — — h): 



N 



AE{Xp,Xh) = eo(Ap)-eo(Ah)+X!(^o(Aj)-eo(Aj)) 



= soiX^) - eo(A^) + f:(£o(Af ) - £o(Af )) + 2^ [ A,^ - Af + f^C^f Xf) 

= Ai?^(Af,Af) + 2^AP^(A^Af), (B7) 
AP(Ap,A;,) = AP^(A^,Af), (B8) 

?B/\B \B\ \ uB(\B \B\ 



where AE {Xp , Xj^) and AP (A^ , A^) are the energy and momentum of a particle-hole excitation in the Bose gas 
with periodic boundary conditions, and A^ and A^ are given by (|B6[) . 

From (|B7[) we see that in the case of twisted boundary conditions, the Fermi velocity uj^*^ will be the same as in 
the Bose gas, whereas for the periodic boundary conditions the Fermi velocity will be modified as 
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